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Residual strains in ferroelectrics are known to adversely affect the material properties by aggravating crack 
growth and fatigue degradation. The primary cause for residual strains is strain incompatibility between 
different microstructural entities. For example, it was shown in polycrystalline ferroelectrics that residual 
strains are caused due to incompatibility between the electric-field-induced strains in grains with different 
crystallographic orientations. However, similar characterization of cause- effect in multidomain 
ferroelectric single crystals is lacking. In this article, we report on the development of plastic residual strains 
in [111] -oriented domain engineered BaTi0 3 single crystals. These internal strains are created due to strain 
incompatibility across 90° domain walls between the differently oriented domains. The average residual 
strains over a large crystal volume measured by in situ neutron diffraction is comparable to previous X-ray 
measurements of localized strains near domain boundaries, but are an order of magnitude lower than 
electric-field-induced residual strains in polycrystalline ferroelectrics. 
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Ferroelectric oxides are used in many electromechanical devices for sensors, actuators, medical imaging and 
energy-harvesting applications, owing to the large values of their piezoelectric coefficients 1 . However several 
reliability issues could be the limiting factors against their applications in many cases. For example, polar- 
ization fatigue and cracking under electromechanical loading have limited applications of ferroelectric single 
crystals in the past 2 ' 3 . Recently, interest in single crystal ferroelectrics has been renewed due to the dramatic 
enhancement of their electromechanical properties with domain-engineering 4 ' 5 . It is therefore imperative that the 
factors that could undermine reliability of domain -engineered ferroelectric crystals be well characterized. It is 
known that in brittle ceramics such as BaTi0 3 , internal crack growth is driven by residual strains generated 
between the different microstructural constituents under the application of electric fields 6 . These internal cracks 
could further aggravate fatigue degradation during repeated electrical cycling 7 . Characterization of residual 
strains in domain-engineered single crystals is therefore desirable. Although generation of residual strains in 
polycrystalline ferroelectric ceramics upon electromechanical loading is well characterized, corresponding stud- 
ies in single crystals has been absent. This is more concerning as large internal strain mismatch is expected 
between the domains of different crystallographic orientations in domain-engineered crystals due to their 
anisotropic piezoelectric properties 8 . In this study, using in situ neutron diffraction we have characterized the 
plastic residual strains and inhomogeneous strain fields in [111] -oriented, domain -engineered BaTi0 3 single 
crystals caused due to interdomain strain incompatibilities during electric field application. 

Domain-engineered ferroelectric crystals are created by the application of a sufficiently large electric field along 
a certain crystallographic axis other than the equilibrium zero -field polar axis. This process creates a set of 
domains whose polarization directions have a common angle with respect to the poling direction 9 . This is 
illustrated in Figure 1(a) for a tetragonal multidomain crystal, where the [111] direction is parallel to the electric 
field and the orthogonal <001> directions of the different possible domain orientations are equally misoriented 
with respect to the applied electric field. Please note that the direction [1 1 1] is different than the equilibrium zero- 
field polar axis of [001] in tetragonal BaTi0 3 . In this sense, the term 'poling' here does not equate to its 
conventional meaning whereby the ferroelectric domains are reoriented so that their [001] axis is parallel to 
the applied electric field; rather it simply refers to the application of a large electric field to a crystal in the as- 
synthesized condition. The exact microstructural changes during this 'poling' process in domain- engineered 
crystals remain unclear, although formation of engineered domain configuration in tetragonal BaTi0 3 has been 
demonstrated for applied electric fields along [111] 10 . Nevertheless, the enhancement in piezoelectric properties 
of domain -engineered crystals has been variously attributed to polarization rotation 11 and/or susceptibility to 
shear deformation 12 . 
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Figure 1 | Crystallographic orientations of domains and experimental set-up: (a) Schematic of domain engineered configuration of [111] -oriented 
single crystal. The pseudocubic [111] direction is equally oriented with respect to the [001] axes of the three orthogonal domain variants, (b) A split-view 
of two adjacent domains across a representative 90° domain wall within [lll]-oriented domain-engineered BaTi0 3 single crystal. Many such domain 
walls can exist within a crystal, a representation of which is shown here. Electric field is applied parallel to 1 1 1 plane normals; while transverse electric- 
field-induced structural changes are indicated parallel to [101] and [110] for Domain 1 and Domain 2, respectively. The misfit in the transverse electric- 
field-induced strain responses of the two domains, as indicated by their d\ x responses, creates an incompatibility in the domain wall region. In the 
schematic, Domain 1 is behind the 90° domain wall and Domain 2 is in front of the 90° domain wall. The spontaneous polarization is parallel to the z 
direction in the two domains, (c) Schematic of in situ neutron diffraction experiment on domain-engineered BaTi0 3 crystal, electric field is applied 
parallel to [111]. All crystallographic coordinates in the schematic are given in pseudocubic notation, (schematic courtesy Renee Manning). 



A particular concern in domain -engineered crystals is the effect of 
strain incompatibilities between the adjacent domains across non- 
180 0 domain boundaries under electromechanical loading. The pre- 
sent work presents a prototypical analysis of this loading scenario, 
wherein a large difference in the electric-field-induced strains exists 
between the two degenerate pseudocubic directions in strongly 
anisotropic BaTi0 3 , as shown in Figure 1(b). In this configuration, 
for both domains, the electric field is applied parallel to the direction 
normal to the (lll)-type planes. For Domain 1, the direction trans- 
verse to the applied electric field is [10l]. Correspondingly, for 
Domain 2, the direction transverse to the applied electric field is 
[110]. In this geometry, the transverse d\ x coefficient corresponding 
to strain parallel to [101] for Domain 1 is —146 pm/V, and the 
transverse d\ x coefficient corresponding to strain parallel to [110] 
for Domain 2 is — 19 pm/V 8 . The notation d\ x denotes the transverse 
coefficient along a certain crystallographic direction such as [110] 
and should be distinguished from the macroscopic transverse piezo- 
electric coefficient of a crystal, d 31 . In domain -engineered BaTi0 3 , 
the large mismatch in the transverse piezoelectric coefficients d\ x of 
the adjacent domains along [110] -type directions is expected to cre- 
ate an electric-field-induced strain mismatch in the domain wall 



region. For an applied electric field of 500 V/mm, the difference in 
d\ x coefficients creates a strain mismatch of —0.006% between the 
domains, which is of the same order of the total electric-field-induced 
strain of the crystal itself. The mismatch between the induced strains 
in domains of two different orientations must be accommodated 
within the microstructure. 

The resolution of strain (in) compatibilities between the different 
components of an aggregate microstructure containing crystallites 
with different orientations is a classic problem in materials science. 
When the adjacent crystallites across a coherent interface exhibit 
strain mismatch as a result of elastic or plastic anisotropics, internal 
stresses are induced to maintain compatibility. Such internally gen- 
erated stresses, which act in addition to external forces, are called 
compatibility stresses 13 . Examples of compatibility stresses can be 
found at interfaces between metallic bicrystals 1314 , at twin boundaries 
in ferroelastic oxides 15 , and at interfaces between epitaxial thin films 
and substrates 16 . When a sufficiently large strain gradient is induced 
across a coherent interface, compensating dislocations are formed 
which locally relieve the compatibility stresses. The compensating 
dislocations, while creating local discontinuities, enhance the overall 
coherency of the interface and minimize the strain energy at the 
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interface due to elastic mismatch 13 " 16 . In ferroelectric materials, 
domain walls are the examples of coherent interfaces wherein com- 
patibility stresses can arise due to difference in piezoelectric res- 
ponses of the adjacent domains 1 . The result of such compatibility 
stresses in domain -engineered crystals is of interest here. 

In the last decade, nanoscale details about the domain structures in 
ferroelectric materials have come forth due to the improved resolu- 
tion of different microscopy techniques. Recent examples include use 
of Piezoresponse Force Microscope (PFM) to reveal the mesoscopic 
dipole closure patterns in single-crystal BaTi0 3 17 , and use of aber- 
ration-corrected transmission electron microscope (TEM) to study 
localized nucleation and growth of ferroelectric domains in BiFe0 3 
thin films 18 . A rather extensive review of recent findings on nanoscale 
domain patterns in ferroelectrics can be found in Catalan et al 19 . 
Although powerful, these microscopy techniques lack quantitative 
information regarding the residual crystallographic strains in ferro- 
electric domain structures, and only indirect information about local 
strain states can be obtained from polarization loops measured with 
PFM 20 . In comparison, scanning diffraction experiments with high- 
energy microfocused or nanofocused X-rays have provided much 
greater insights regarding the strain fields around domain walls in 
ferroelectric single crystals 21 " 23 . However, scanning based techniques 
with either pointed tips or X-rays are not suitable for dynamic mea- 
surements during in situ electric field application because the posi- 
tion of interest changes with domain wall motion 24 . Moreover, 
scanning diffraction experiments can only provide strain informa- 
tion at a particular local point and elucidating residual strains over 
larger sample volumes would require a statistical analysis of several 
different measurements. Therefore, it would be much more advant- 
ageous if internal strains in ferroelectric domain structures could be 
revealed from diffraction measurements taken with a large beam that 
is comparable in size with the material of interest. 

In this study our objective is to investigate the evolution of internal 
strain fields in [1 1 1] -oriented and domain -engineered BaTi0 3 ferro- 
electric crystals during and subsequent to the application of electric 
fields. The strain fields arise due to alterations in the atomic-scale 
structure that are readily measured in diffraction. In situ neutron 
diffraction of single crystals under applied electric fields provides 
volume-averaged information of strain states across multiple 
domains and the evolution of internal strain fields can be ascertained 
from such information, as shown here. The particular advantage of 
neutron diffraction is that neutrons penetrate deep inside a sample 
and therefore probes structural changes inside the bulk of a material. 
BaTi0 3 is chosen in the present work to serve as a prototypical 
ferroelectric in which to observe the effects of strain mismatch 
and, moreover, to avoid the structural complexities that could arise 
in relaxor ferroelectrics due to a coexistence of polar nanoregions and 
multiple ferroelectric phases 25 " 27 . 

Results 

The experimental setup for in situ diffraction under applied electric 
fields is illustrated in Figure 1(c). The BaTi0 3 crystals were obtained 
in their as-processed state and subjected to the electric fields for the 
first time during the current experiment. Two different crystals were 
used, one of which was subjected to positive electric fields and 
another to negative electric fields. The neutrons diffracted from the 
sample were measured on 2-D area detectors (see Supplementary 
Information for details). All in situ neutron diffraction measure- 
ments were taken for applied field amplitudes of greater than (+/ 
— ) 100 V/mm, which is the nominal coercive field for tetragonal 
BaTi0 3 single crystals 10 ' 28 . 

Structural changes measured parallel to applied electric field. 

Representative changes in the integrated 333 diffraction peak 
profile upon application of positive and negative electric fields are 
shown in Figures 2(a) and 2(b), respectively. The 333 momentum 



transfer vector for neutron diffraction is parallel to the direction of 
the applied electric field, and therefore provides measurement of 
crystallographic changes parallel to the applied field direction. In 
both cases, with increase in field amplitude, the 333 peak position 
shifts towards higher <i-spacing, along with an increase in total peak 
intensity. Electric-field-induced changes in the (333) interplanar d- 
spacing, Ad, are calculated with respect to the d-spacing measured at 
a certain minimum field amplitude of ( + /— ) 100 V/mm. As shown 
in Figure 2(c), Ad increases linearly with increase in electric field 
amplitude for both samples, except for the field amplitude of 
275 V/mm, at which point a sharp discontinuity in Ad is evident 
for the sample subjected to positive electric fields. No such 
discontinuity is observed for the crystal subjected to negative 
electric fields. Linear fits to Ad yield corresponding electric-field- 
induced strain coefficients of 678 pm/V and 845 pm/V, for the 
crystals subjected to positive and negative electric fields, 
respectively. These values of strain coefficients are more than 3 
times the values of d^ 3 parallel to the 111 crystallographic direction 
(—200 pm/V), as calculated from the coefficients of single crystalline 
piezoelectric tensor of BaTi0 3 12 or reported from macroscopic 
measurements 29 . This indicates that the measured electric-field- 
induced 333 lattice strains observed here are more complex than 
the intrinsic piezoelectric strain of the lattice suggests. 

The normalized 333 peak intensities are plotted as Af u /A° lu in 
Figure 2(d), where Af n refers to the integrated 333 peak intensity for 
an applied electric field amplitude E and A° 1U refers to the integrated 
333 peak intensity for an applied electric field amplitude of (+/ — ) 
100 V/mm. Figure 2(d) shows that the 333 peak intensity increases 
linearly with an increase in the amplitude of negatively polarized 
electric field. However, the crystal subjected to positive electric field 
shows a minimal change by comparison. Notably, a sharp discon- 
tinuity in the 333 integrated peak intensity can be observed at the 
field amplitude of 275 V/mm, which correlates with the trend 
observed for Ad for the crystal subjected to positive electric fields. 

Structural changes measured transverse to applied electric field. 

Representative changes in 220 diffraction peak profiles upon 
application of electric fields are shown in Figures 3(a) and 3(b). 
The 220 momentum transfer vector lies orthogonal to the 
direction of the applied electric field, and therefore provides 
measurement of crystallographic changes transverse to the applied 
field direction. The results for the normalized peak intensities are 
represented as A^q/A^q in Figure 3(c), where the superscripts have 
the same meaning as for the normalized 333 diffraction peak 
intensities. The change in integrated 220 peak intensities is larger 
for the crystal subjected to negative electric fields as compared to the 
crystal subjected to positive electric fields. Also, consistent with 
the trend observed in Figure 2(d), Figure 3(c) shows a sharp 
discontinuity in the integrated 220 peak intensity at 275 V/mm for 
the sample subjected to positive electric field. It is important to note 
that the observed intensity changes for the 220 diffraction peaks are 
different from intensity interchange observed normally between 
degenerate diffraction peaks in ferroelectric materials that are 
characteristic of change in volume fraction of non-180 0 domains 30 . 
In fact, the peak corresponding to the reflecting planes normal to the 
[220] direction exhibits greater change in intensity as compared to 
the peak corresponding to the reflecting planes normal to the [202] 
direction. Therefore, it is unlikely that the intensity change for the 
220 diffraction peaks observed here is due to non-180 0 domain 
switching. The origins of the measured lattice strains and peak 
intensity changes are examined in the following discussion. 

Discussion 

As shown in the schematic in Figure 1(b), the direction normal to the 
(hhh) planes is parallel to the (kkO) -type planes that correspond to 
the 90° domain walls 29 . Therefore, during the application of electric 
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Figure 2 | Electric-field-induced structural changes parallel to direction of applied electric field: (a),(b) changes in 333 diffraction peak profiles for the 
crystal subjected to (a) positive, and (b) negative electric fields; (c) change in (333) lattice plane spacing, Ad, as function of applied electric field 
amplitude; (d) change in integrated intensity of 333 diffraction peaks as function of electric field amplitude. The arrows indicate the directions of 
progressive changes in the applied electric field. The final datapoint is measured at zero field subsequent to the application of a large electric field. 



fields, the internal compatibility stresses arising from a strain mis- 
match across 90° domain walls causes straining of the (hhh)-type 
planes, which is evident from Figure 3(c). The corresponding lattice 
strains are calculated from the data shown in Figure 3. A maximum 
(hhh) lattice strain (Ad/d) of ~2 X 10~ 4 is observed for an applied 
field amplitude of 400 V/mm. Additionally, subsequent to the elec- 
tric field application, a residual strain of ~7 X 10" 5 is observed for 
the crystal subjected to positive electric fields, and a corresponding 
residual strain of ~2 X 10" 4 is observed for the crystal subjected to 
negative electric fields. These values of residual (hhh) lattice strains 
are comparable to the zero-field (002) lattice strains of ~7 X 10" 5 
measured with microfocused X-rays around 90° domain walls in 
tetragonal BaTi0 3 crystals 21 . In comparison, residual (hhh) lattice 
strains of ~10~ 3 were reported for polycrystalline lead zirconate 
titanate (PZT) ceramics subsequent to electric field application 31 . A 
higher residual strain in polycrystalline PZT could be due to the 
higher piezoelectric coefficients of soft PZT materials as compared 
to that of tetragonal BaTi0 3 . On the other hand, while the piezo- 
electric coefficients of PZT are higher by about a factor of five, the 
residual lattice strains measured in polycrystalline PZT is more than 
an order of magnitude higher than the corresponding lattice strain 
values measured here for BaTi0 3 single crystals. This likely suggests 
that the presence of grain boundaries enhances the residual strains in 
ferroelectric materials. 

It was observed earlier in polycrystalline tetragonal ferroelectrics 
that strain mismatch between the neighboring grains causes residual 
elastic strains for the (hhh) -type lattice planes 31 . The residual hhh 
elastic lattice strain for a set of similarly oriented grains (measured 



subsequent to electric field application) can be either positive or 
negative, depending on their specific crystallographic orientations. 
However, the net elastic strain integrated over all the grains vanishes 
due to the stress balance condition and the assumed uniform elastic 
constant 32 . This situation can be compared with the present case 
wherein a strain mismatch exists between the differently oriented 
domains in a single crystal sample. For single crystal diffraction, 
the hhh poles from all the domains are measured along a single 
direction. Therefore, the overall measured 333 lattice strain is a 
volume-average of the 333 lattice strains of all the domains, which 
includes domains on both sides of the 90° domain walls. In this 
situation, based on internal stress balance condition, no residual 
elastic lattice strains should be observed at zero electric field if the 
internal strain fields across 90° domain walls are symmetric. The fact 
that a residual hhh elastic lattice strain is measured subsequent to 
electric field application indicates a generation of asymmetric 
internal strain fields across 90° domain walls. Moreover, it also shows 
a certain degree of plastic deformation within the crystal. 

In addition to changes in lattice d- spacing, internal gradients of 
lattice distortions in single crystal samples also result in an increase of 
diffraction peak intensities, as was previously observed in ferroelec- 
tric crystals upon application of electric fields 27 ' 33 . The currently 
observed changes in diffraction peak intensities can be similarly 
attributed to the evolution of internal strain gradients. The differ- 
ently oriented domains with orthogonal polarization vectors across 
90° domain walls give rise to the two degenerate Bragg peaks for the 
220 reflection. It is observed here that as an effect of electric field 
application, these two peaks exhibit disproportionate changes in 
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Figure 3 | Electric-field-induced structural changes transverse to the 
direction of applied electric field: changes in 220 diffraction peak profile 
for crystals subjected to (a) positive, and (b) negative electric fields; (c) 
change in integrated intensity of 220 diffraction peaks as function of 
electric field amplitude. The arrows indicate the directions of progressive 
changes in the applied electric field. The final datapoint is measured at zero 
field subsequent to the application of a large electric field. 

peak intensities (Figure 3(a) and (b)). This indicates that the internal 
strain gradients on either side of the 90° domain walls are not sym- 
metric. This view is consistent with the earlier assertion that asym- 
metric strain fields across 90° domain walls are also the cause for the 
observed residual 333 lattice strains. In fact, we observe an exact 
correlation between the residual 333 lattice strains and the increase 
in diffraction peak intensities - the sample that shows a larger resid- 
ual strain also exhibits a larger increase in 220 and 333 peak intens- 
ities. A hypothetical case where a 90° domain wall can have an 
asymmetric strain profile is schematically illustrated in Figure 4. 

Previous microdiffraction experiments on single crystals indeed 
showed that some components of deviatoric strain tensors exhibit 



asymmetry as a function of position around a 90° domain wall, even 
under zero electric fields 22 . Such strain asymmetries were also 
reported along twin boundaries in germanium bicrystals due to the 
presence of dislocations 34 . Strain asymmetries can arise due to a 
difference in the stiffness of the compressed and expanded regions 
around a dislocation along the twin boundaries, which in turn 
depends on the local atomic arrangements at the interface 34 . It can 
be debated whether dislocations are also responsible for asymmetric 
strain fields across 90° domain walls, although that scenario seems 
likely in the present case. Since the transverse strains of the adjacent 
domains differ by a factor of ~7, there is a large strain mismatch in 
the 90° domain wall area and this consequently increases the local 
strain energy. The strain energy near 90° domain walls can be poten- 
tially reduced due to interaction between the strain fields of disloca- 
tions with Burger's vector bll (110) and those of 90° domain walls, as 
recently shown by phase-field modeling 35 . It is therefore possible that 
the increased strain energy due to interdomain strain incompatibility 
is partially relieved by new dislocations created along 90° domain 
walls. This is one possible explanation for the electric-field-induced 
residual strains in domain -engineered BaTi0 3 crystals. Additionally, 
the inhomogeneous field distribution near the 90 °C domain wall due 
to non-unifom dielectric response could also be a contributing factor 
for strain asymmetries near 90 °C domain walls. For further insights, 
it will be interesting to look at the microscopic length scales over 
which the internal strains are distributed around the domain wall 
region in domain-engineered crystals. It is expected that future 
experiments with micro/nano focused X-rays can provide such 
information. 

The reason for the slightly different strain behavior of the two 
samples is not entirely clear, although multiple effects could be con- 
tributing to this such as extension (sample subjected to negative 
electric field) vs. compression (sample subjected to positive electric 
field) of the crystal in the transverse direction, or presence of pre- 
existing defects along the domain walls. Nevertheless, the subtle 
differences in the structural changes observed for the two different 
samples suggest that the electric-field-induced structural changes in 
ferroelectric single crystals are highly sensitive to the existent local 
imperfections and their evolution under applied electric fields. 

The connection between interdomain residual strains and per- 
formance of ferroelectric single crystals is an area that needs further 




Twin boundary 



Figure 4 | Domain-wall with asymmetric strain profile: The proposed 
structure of a 90 domain wall with an asymmetric strain profile across it. 
The arrows indicate the direction of polarization vector for the two 
domains. 
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clarification. Future research efforts should be directed to under- 
stand effects of residual strains on time -dependent and fatigue beha- 
vior of multidomain single-crystalline functional materials. 

Methods 

Two single crystals of BaTi0 3 in the engineered domain configuration were obtained 
from Mateck Gmbh. The two crystals were obtained from two different processing 
batches. The crystals were cut so that their face normals were parallel to the ortho- 
gonal (111), (lOl) and (121) directions, and electric field was applied across the 
crystal faces normal to the (111) direction. In situ structural measurements during the 
application of electric fields were obtained at the newly commissioned VULCAN 
Engineering Diffractometer of the Spallation Neutron Source at the Oak Ridge 
National Laboratory 36 . The structural changes in the two different crystals under 
opposite polarities of applied electric fields were measured. In one sample, the 
amplitude of the positively polarized electric field was increased stepwise with 
increments of 25 V/mm, starting from ( + ) 100 V/mm up to a maximum field of ( + ) 
400 V/mm. In the second sample, the amplitude of the negatively polarized electric 
field was increased from ( — ) 100 V/mm up to a maximum field of ( — ) 400 V/mm, in 
steps of 25 V/mm. Neutrons diffracted from the crystal are measured by the two area 
detectors placed at 1 80° degrees with respect to each other. The vertical axes of the two 
detectors in the laboratory 3-D space are parallel to the (112) crystallographic 
directions, while the horizontal axes of the detectors are in the plane containing the 
(lll)and (101) the crystallographic directions. Therefore, the two area detectors of 
the diffractometer measure diffraction peaks in the time-of-flight mode, with their 
momentum transfer vectors approximately parallel to the orthogonal (ill) and (101) 
directions, respectively. The measured intensities of the diffracted neutrons are 
summed over multiple pixels that have the same coordinate along the horizontal 
detector axis. For details of the measured diffraction intensities, please see 
Supplemental information. 
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